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Abstract 

A semiclassical theory of dissipative Henon-Heiles system is proposed. Based 
on ^-scaling of an equation for evolution of Wigner quasiprobability distri- 
bution function in presence of dissipation and thermal diffusion, we derive a 
semiclassical equation for quantum fluctuations, governed by the dissipation 
and the curvature of the classical potential. We show how the initial quan- 
tum noise gets amplified by classical chaotic diffusion which is expressible in 
terms of correlation of stochastic fluctuations of the curvature of the potential 
due to classical chaos and ultimately settles down to equilibrium under the 
influence of dissipation. We also establish that there exists a critical limit to 
the expansion of phase space. The limit is set by chaotic diffusion and dissi- 
pation. Our semiclassical analysis is corroborated by numerical simulation of 
quantum operator master equation. 
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I. INTRODUCTION 



The influence of dissipation on quantum dynamics of classically chaotic systems has 
been one of the key issues in nonlinear dynamics today. The dynamical way of dealing 
with dissipation is to consider a system-heat bath model which has been the cornerstone 
for understanding dissipative processes* 1 ) in a wide range of physical disciplines^ 2 ) , such as, 
condensed matter physics, quantum optics, chemical dynamics etc. The theoretical develop- 
ment in this regard is welldocumented in the literature*- 2-6 -*. When the system, in question, 
is classically chaotic, one envisages a variety of rich physics^ 7-20 ) concerning localization 
and its suppression, quantum measurement problem, irreversibility, relaxation, decoherence 
etc. The inferences drawn from these are sometimes extended to the question of generic 
quantum chaos. For example, a dissipative quantum system exhibting chaos in its classical 
limit was constructed by coupling the quantum kicked rotor to a reservoir to obtain dissi- 
pative quantum standard map by Dittrich and Graham ( - 7 \ It was observed that even weak 
damping is capable of disrupting dynamical localization which suppresses chaotic motion in 
the conservative standard map and thus restores diffusion in action variable in the timescale 
of classical relaxation. The effect of quantum correlation on classical chaotic behaviour 
had been illustrated by Sundaram and Milonni® by considering a kicked quantum system 
coupled to a reservoir. An appropriate choice of potential results in a logistic map with 
self-consistently generated quantum correlations. It has been observed that at intermediate 
range of dissipation an irregular behaviour is induced by quantum correlations even when the 
classical limit is regular. Based on an analysis of quantum Brownian motion in d-dimensions 
using the unified model for diffusion localization and dissipation, Choem 20 ) has proposed a 
semiclassical strategy at low temperature using Feynman- Vernon propagator scheme. It has 
been demonstrated that different mechanisms for dephasing emerge for ergodic and noner- 
godic motions. In another issue Bonilla and Guinea^ 15 -* have studied a simple model having 
quantum and classical degrees of freedom in presence of dissipation. The emergence of chaos 
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in an open quantum system has also been considered by Spiller and Ralph* 16 ). 

The majority of the studies considered above are based on maps, (such as, standard 
map or logistic map) which have been the testing ground for various theories of chaos. 
We construct here a dissipative version of a two-degree-of-freedom continuous system - the 
Henon-Heiles model ( 21_23 \ to study the evolution of a quantum system in presence of dissi- 
pation and thermal diffusion. The Henon-Heiles model captures the essential generic features 
of classical chaos in nonintegrable systems and has been widely applied in the context of 
astronomy and chemical dynamics over the last several decades^ 21-23 ). Based on suitable 
^-scaling of Wigner equation which incorporates the effect of dissipation and thermal diffu- 
sion, we formulate a semiclassical dynamics which is governed by dissipation and curvature 
of the classical potential. The stability of classical motion is determined by the nature of 
curvature of the potential which in the chaotic regime can be considered to be a stochastic 
process^ 24 ' 25 ) . An appropriate treatment of this stochastic process in terms of the theory 
of multiplicative noise yields a Fokker-Planck equation of motion for Wigner-function. We 
design the initial conditions in terms of minimum uncertainty wave packets to maximize the 
classical-quantum correspondence and show how the initial wave packet corresponding to a 
chaotic trajectory evolves in time, and how the initial quantum noise (inherent in minimum 
uncertainty of wave packet) associated with it gets amplified by intrinsic classical stochastic- 
ity at the begining to eventually settle down to equilibrium under the influence of dissipation. 
We establish that there exists a critical limit to the expansion of phase space. The limit is set 
by chaotic diffusion and dissipation. Our semiclassical analysis is supplemented by quantum 
simulation of the operator master equation to verify the basic theoretical propositions. 

The outlay of the paper is as follow: In Sec. II we introduce the quantum operator 
master equation and the Wigner function equation for an open system, ^-scaling of the 
Wigner equation results in a semiclassical equation governed by the dissipation due to the 
surroundings and the curvature of the potential. This is followed by van Kampen's treatment 
of multiplicative noise ( 26 ' 27 ) to deal with stochastic fluctuations of the curvature of the 
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potential which leads to a Fokker-Planck equation. In Sec III the Fokker-Planck equation 
is adapted to Henon-Heiles system followed by a detailed analysis of the problem. In Sec 
IV a numerical simulation of the operator master equation has been carried out to verify 
the theorerical propositions. The approximations and their validity with a summary of the 
main results have been discussed in Sec V. 



II. CHAOTIC EVOLUTION OF AN OPEN SYSTEM; GENERAL ASPECTS 



A. Quantum dynamics 



To study^ 2 ) the evolution of a quantum system in presence of weak dissipation and 
thermal diffusion we first consider the Hamiltonian of an N-degree-of-freedom system H . 

N 2 

Ho = J2£r + V(tii}) , i = l--N (1) 

i=l 1 

where {qi,Pi} represents the coordinates and momenta of the N-degree-of-freedom system. 
V({qi}) is a nonlinear potential such that the classical version of H admits of chaos. 

The bare system is then coupled to an environment modeled by a reservoir of harmonic 
oscillators, governed by the following total Hamiltonian 

oo oo 

H = H + h^2 tijbj% + hY, [k(£lj)bj + k*{Qj)b^] q , (2) 

3 3 

where b 3 (bf) denotes the annihilation (creation) operator for the harmonic oscillators 
which comprise a bath. The third term represents the linear coupling of one of the selected 
degrees of freedom (through co-ordinate q) of the system to the bath. k(Qj) signifies the 
coupling constant. 

It is convenient to invoke first the rotating wave approximation (RWA). After appropriate 
elimination of reservoir variables in the usual way using Born and Markov approximations 
we are lead to the following standard reduced density matrix equation for the evolution of 

the system^ 2 ), 



— = — — [Ho, p] + — (2apa' — a'ap — pa) a) + D(a) pa + apa' — a'ap — paa') . (3) 

Here we have expressed the system operators q,p (for N — 1) in terms of a harmonic 
oscillator operators a(annihilation) and a' as g = :^=( a + ^) an d p = i^J !1 f i (a — a)). Note 
that the harmonic oscillator characterized by frequency u has nothing to do with the reservoir 
of harmonic oscillators. In the derivation above, one uses a broad band spectral density 
function for the reservoir evaluated at uj to realize the damping constant 7 as 2T\-\k{oj)\ 2 g{ui) 
within a Markovian scheme. D(= fry) is the diffusion coefficient and n(= [exp (jjfj — 
refers to the average thermal photon number of the reservoir. The terms analogous to Stark 
and Lamb shifts are neglected. If more than one degree of freedom of the system is coupled 
to the bath then the coupling term in Eq. (2) and dissipative terms (7 and D terms) in Eq. (3) 
should appropriately include additional similar contributions [see Eq.(14)]. 

The first term in Eq.(3) corresponds to the dynamical motion of the system that generates 
Liouville flow and the second term denotes the loss of energy from the system to the reservoir, 
while the last term indicates the diffusion of fluctuations of the reservoir into the system 
of interest. The terms containing 7 arise due to the interaction of the system with the 
surroundings. 

We note that Eq.(3) is a popular form of the operator master equation (as derived by 
Louisell^ 2 )) which is widely used in quantum optics. The equation has also been applied 
earlier by Dittrich and Graham^ 7 ) in the treatment of dissipative standard map and the 
related problems of chaotic dynamics by others'^ 8 ' 16 ). The Eq.(3) is also general in the sense 
that we need not ascribe any notion of regularity or chaoticity in describing the motion gov- 
erned by the Hamiltonian system(iJo)- The correlation between different forms of operator 
master equation has been reviewed in Ref^ 5 ^. All of them, however, are not well-suited for 
numerical simulation. Eq.(3) suits this purpose well. We shall return to this issue in Sec. IV 
to verify the theoretical propositions. We note, in passing, that Eq.(3) is based on rotating 
wave approximation and Born-Markov approximation. The latter approximation restricts 
its validity to weak damping limit only. 
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B. Semiclassical theory 



Our next task is to go over from a full quantum operator problem to an equivalent c- 
number problem described by the Hamiltonian(2). To this end we consider the quasi-classical 
distribution function W ({(&}, {Pi}, t) of Wigner^ 28 ). The time evolution of this phase space 
function of the dynamical system characterized by the c-number variables {qi,Pi} is based 
on two considerations: First, one takes into account of the usual dynamical evolution under 
the influence of potential V as defined in (1). The second is the dissipative evolution of the 
system when it is coupled to the harmonic oscillator bath described by Eq.(2). The former 
is essentially rewriting Schrodinger equation in a quasi-classical langauge and has nothing 
to do with the latter. Thus we write 

(dW\ _ (dW\ (dW\ 
While the dynamical evolution is governed by Wigner equation^ 28 ) , 



'dW\ 



N 



dt ) dynamcal S 



Pi dW , (dV\ dW 



2rrii dqi \ % J dpi 



+ J2 ( &*+- n »v \ (!) 



mH hnjv-1 



dqT---dq n N N ) n 1 \---n N \ 



n i + n 3 + ■ ■ ■ + n N is odd and > 1 
gm+---+N N 



x W 

dpT---dp n N N ' 

the form of ( ^ ) is due to Caldeira and Leggett^ as given by [when one of the 

V I dissipative 

system degrees of freedom is coupled to the reservoir as expressed in Hamiltonian (2); see 
the dissipative part of Eq.5.14 of Ref(5)] 



/ dW\ n d d 2 W 

= 2^—pW + D-—- 

dissipative 



dt j ,. . ,. " 1 drF dp 2 

/ dissipative r r 

where 7 and D have the same significance as in Eq.(3). The first term in the last equation 
is a direct consequence of the existence of a 7-dependent term in the imaginary part of the 
exponent in the expression for the propagator for the density operator of Feynman and 
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Vernon theory and has been shown^ 5 ) to be responsible for appearance of a damping force 
in the classical equation of motion for the Brownian particle to ensure quantum-classical 
correspondence. 

The total dynamics is a superposion of two contributions provided by the last two equa- 
tions and when written elaborately we have; 



dW_" 



Pi dW (dV\ dW 



2rrii dqi \ dqi I dpi 



( d n 1+ -n N y 

+ E 



niH hriAr — 1 



dqT---dq n N N ) n 1 \---n N \ 



n>i + ris + ■ ■ ■ + tin is odd and > 1 

d n 1+ -+N N Q Q2 W 

W + 27 — vW + D 

dpT---dp7 dp dp 2 



x ,^ N W + 2^-pW + D-—- . (4) 



That the two contributions to the total evolution of the Wigner function in Eq.(4) act 
independently in the overall dynamics is an assumption. This assumption is also implicit 
in the operator master equation (3) and has been routinely used in nonlinear and quantum 
optics, in general. (Note that Eq.(5.14) of Ref(5) carry the same messege for a single-degree- 
of-freedom system). Strictly speaking, the 7 and D terms in Eqs.(3) and (4) are valid if 
the system operators [i. e. , q and p in Eq.(2)] pertain to a harmonic oscillator. When the 
system is nonlinear, as in the present case (also in many nonlinear optical situations) the 
usual practice is to add the additional contribution — i[H non , p] to the master equation [ in 
the language of Fokker-Planck description this commutator, in general, contributes higher 
(third or more) order derivatives of the distribution] and to assume that the dissipative 
terms remain unaffected by the addition of the commutator term, H non being the nonlinear 
part of the Hamiltonian H . The validity of this assumption was examined^ 29 ) earlier by 
Haake et. al. and also by us. It is now known that this assumption is quite satisfactory 
within the perview of weak damping and/ high temperature limit. 

The equation (4) is a full quantum mechanical equation. A simple version of the above 
equation for one-degree-of-freedom system was used earlier by Zurek and Paz^ 18 ) for studying 
some interesting aspects of quantum-classical correspondence in relation to decoherence. 



The primary reasons for choosing Eq.(4) as our starting point for semiclassical analysis 
are (i) the rotating wave approximation (RWA) in the system-reservoir coupling has not 
been made in deriving Eq.(4) (ii) Eq.(4) is also free from Born approximation (or weak- 
coupling approximation) ensuring that the theory is valid even in the strong damping limit, 
(hi) Eq.(4) reaches the correct classical limit when h — > 0, when D becomes the thermal 
diffusion coefficient in the high temperature limit. Thus Eq.(4) is a good description in the 
semiclassical limit. Keeping in view of these remarks and the earlier discussion in Sec IIA we 
observe that Eq.(3) and Eq.(4) describe same dynamics in weakly dissipative systems. We 
adopt Eq.(4) for our semiclassical analysis that follows and Eq.(3) for quantum numerical 
simulation in Sec IV to verify the theoretical propositions of this analysis. 

The first term in Eq.(4) is the usual Poisson bracket which generates the Liouville flow. 
Both the Poisson bracket and the higher derivative terms result from an expansion of the 
Moyal bracket on the basis of an analytic V(q). The last two terms have the same significance 
as in Eq.(3). It is important to note that the failure of correspondence between classical 
and quantum dynamics is predominantly due to higher derivative terms which make their 
presence felt roughly beyond the Ehrenfest regime. 

As a first step it is convenient to introduce the following scaling of c-numbers {qi,Pi} in 
analogy to van Kampen's f2-expansion; 

ft = Qi(t) + h 1/2 /j,i , 

Pi=Pi{t)+h 1/2 Vi , (5) 

where h is the associated smallness parameter for the present analysis, /i and v in 
Eq.(5) refer to quantum fluctuations in co-ordinate and momentum, respectively. q(t) and 
p{t) are the corresponding classical co-ordinate and momentum. The time evolution of the 
distribution function of the fluctuation variables obeys 



dt " " V 



v k dcj) d 2 V dcj) 



+ 2 7 |-W> + Oih 1 ' 2 ) . (6) 



m k d/j, k dqjdq k dis k 
Although this equation does not involve any h explicitly, it describes the time evolution 

of probability density function 0({/i«} , {^} , t) for the quantum noise variables {fa, i/j}, since 



is the lowest order quantum correction to classical distribution function W (qi(t),Pi(t),t). 
Secondly, the quantum dynamics enters into the picture when we put the quantum con- 
straint^) on the initial density function (fi({fik} , {^k} , 0) 



N 



<K{M , {vk\ , o) = n t~ ex p 
k=i 4<t 



0^2,, 2 



(7) 



as 



1 



1 



(8) 



where h — 1 is used. 

We thus note that the initial density z/, 0) is not a 5-function but has an appropriate 
spread. This spread incorporates the quantum noise which gets amplified as the density 
evolves in time. It is thus a quantum (minimum uncertainty product) condition and a 
requirement imposed by quantum-classical correspondence, v in Eq. 6 refers to the specific 
degrre of freedom of the system to which the reservoir is coupled to allow the exchange of 
energy between the system and the reservoir. 

As a second step we put Eq.(6) in a more compact form by invoking the symplectic 
structure of the Hamiltonian dynamics. For this, we specify 



qi for % — 1 • • • TV , 
Pi- N fori = N + l,---2N . 



(9) 



Defining I as 



I = 



( o e\ 



\ 



—E 

where E is an iV ® N unit matrix, one can write the Hamilton's equation 



(10) 



(11) 



Again we introduce the scaling 



U = Zi(t) +h 1/2 rn 



(12) 



with 

rji = ^ for i = 1, • • -N 

= Vi_ N for i = N + 1 • • • 2N , (13) 

corresponding to quantum fluctuations in co-ordinates (/!«) and momenta (i/j). We gener- 
alize Eq.(6) further to the extent that all the momentum components (rji, i = N+l, ■ ■ ■ 2N) 
are coupled to the bath linearly. One obtains the equation of motion for quantum fluctuation 
distribution function 

d 



<90 v 

hJ 



(14) 



where we have assumed that 

li = for j = l,---N , 

7,- = 7 for j = + 1, • • -2N . 

Here 

J - = ^ fc ^ (15) 
contains the second derivatives of the potential and is a function of classical dynamical 
variables Zi(t), (i. e., Pi(t) and <&(£)). 

For further treatment Eq. (14) may be rewritten in a more compact form as follows: 

^ = [-F(f) • V + 27V 7 ]0 , (16) 

where 

F(t) = J(t)r] - 2-fKri . (17) 

V refers to differentiation with respect to the components of rj and K is a 2N <g> 2N matrix 
defined as 

kij = for % ^ j , 

ku — for « = 1, • • • TV , 

A- = 1 for i = N + 1,---2N . 
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J, the jacobian matrix as defined in (15) is a function of classical dynamical variables 
{qi(t) , Pi(t)} . The crucial question of stability/instability of classical motion in Hamiltonian 
systems essentially rests on this jacobian, or curvature (or second derivative) of the poten- 
tial. Traditionally the local linear stability analysis around the fixed points is based on the 
assumption^ 22,23 ) of constant curvature. However, the true stability of motion is only deter- 
mined by keeping the time dependence of J (implicitly through , Pi(t)}) matrix intact. 
Also there is little connection between the local stability and global chaos. In view of this 
it is necessary to take full account of the time dependence of the curvature of the potential 
J(t) along the trajectory itself. When the motion of the dynamical system is regular J(t) is 
highly correlated throughout the entire course of evolution. On the other hand for chaotic 
motion when the dynamical variables in J (i. e. , {qi(t) , Pi(t)}) by virtue of the classical 
equation of motions for q^t) and Pi(t) [or in general Zi(t) of Eq.(9)] behaves stochastically 
J(t) describes a stochastic process. The loss of correlation in chaotic dynamical systems 
thus rests on the decay of correlation of fluctuations of J(t). What follows subsequently is 
a stochastic description of classical chaos in terms of this correlation. 

Ever since the early numerical study of Chirikov mapping^ 30 ) revealed that the motion 
of a phase space variable {q or p} can be characterized by a simple random walk diffusion 
equation, attempts have been made to describe chaos in terms of a stochastic description ( 
Langevin and Fokker-Planck description has been widely employed). It has now been realized 
that deterministic maps can result in long time diffusional processes and methods have been 
developed to predict successfully the corresponding diffusion coefficients'- 31 ^ . In a number of 
recent studies^ 11 ' 24,25 ) we have shown that the fluctuation in the curvature of the potential 
is amenable to a stochastic description in terms of the theory of multiplicative noise. This 
allows us to realize a number of important results of nonequilibrium statistical mechanics, 
like Kubo relation'- 24 -' fluctuation-decoherence relation'- 25 -' etc. in chaotic dynamics of a few- 
degree-of-freedom system. 

Another important point to be noted here is that we do not make any a priori assumption 
about the nature of the stochastic process (J(t)). The special cases, such as, noise is Gaussian 
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or 5-correlated, etc. have attracted so much attention in the literature that it is necessary to 
emphasize that we have not made any such approximation. The stochasticity of F(t) depends 
on J which is determined by the exact solution of the classical equation of motion(ll). 
Eq.(16) may therefore be regarded as a stochastic differential equation with multiplicative 
noise. For convenience F(£) • V can be partitioned (this partitioning will be clarified in more 
detail in the next section) into two parts; a constant part F ■ V and a fluctuating part 
Fi(f) • V. Thus we write 

F V = F - V + Fi-V . (18) 

We now come to the third step. Making use of one of the main results for the theory of 
linear equation of the form (16) with multiplicative noise, we derive an average equation for 
as given by (for details, we refer W 26 ' 27 )); 

5(0) 



dt 



= j-F • V + 2iV 7 - (Fx • V) + dr«Fi(i) • V exp(-r[F • V + 2iV 7 ]) 

Fl (t - r) • V» exp(r[F • V + 2iV 7 ])} (0) (19) 



where ((• • •)) implies ((<Mj)) = (<Mj) — (?.?)■ The operator exp(— rF • V) provides 
the solution of the equation 

^M = -F .V/(^0 (20) 

(where / signifies the "unperturbed" part of 0) which can be found explicitly in terms 
of characteristic curves. The equation 

rj = F (77) (21) 

for fixed t determines a unperturbed mapping from i](t = 0) to t)(t), i. e., 77 — > 77^ with 
inverse (rj T )~ T = 77, The solution of (20) is 

drj' 



= expHFo-V]/M). (22) 



drj 1 
drj 



drj 

being a Jacobian determinant. The effect of exp(— tF • V) on f(rf) is as follows 
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c M-tF -v)f( v ,o) = f( v - t ,o) 

When this simplification is used in Eq.19 we obtain 

<9(0> 



drf 



t , - | -F .V + 2iV 7 - (Fx. V) + J™ dr 
((FxfoO.VrFxfa-V-r^-V-, 



drj 
drf T 



(23) 



d?7 
drj 



drf 



(0) • 



(24) 



The above consideration is based on a second order expansion in ar c (by van Kampen^ 26 ) ) , 
where a is the strength parameter required for bookkeeping the order of the perturbation 
fluctuation and r c is the correlation time of fluctuations in Fi(t) [in the derivation above 
we have put a — 1]. The average (0) in Eq.(24) varies on a coarse-grained timescale 
which is much slower compared to the timescale set by the correlation time of fluctuation of 
Fx(t). Second, the derivation above neglects the effects of higher powers of U and thus the 
Eq.(24) is an effective semiclassical equation for quantum fluctuation distribution function. 
Since it contains second derivatives with respect to components of 77, it has the form of a 
Fokker-Planck equation. Third, the theory discussed so far (Eq.(24)) is valid, in general, for 
N-degree-of-freedom systems. 

We now adapt Eq.(24) to the classic paradigm of chaotic dynamics - the Henon-Heiles 
system. 

III. THE DISSIPATIVE HENON-HEILES SYSTEM 

A. The Fokker-Planck equation 

We consider the Henon-Heiles system which is kept in contact with the surroundings. 
The Hamiltanian of this system is given by 



Ho = £ + f + V( qi ,q 2 ) , 



(25) 



where V(qi, q 2 ) = \{q\ + q\ + 2qfg 2 — |gf ), is the potential energy of the two-degree-of- 
freedom system. 
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The classical equations of motion of the particle in presence of damping (at a rate 7) are 



Qi = Pi , 

dV{q 1 ,q 2 ) . 
Pi = -lPi ~ ,? = 1,2 



(26) 



Note that in the above equation we have assumed for simplicity the value of dissipation 
rate same for both the degrees of freedom. The equations of motion for the quantum 
fluctuation variables 771 ,772, i] 3 and 774 corresponding to qi, q 2 , Pi and p 2 , respectively, read 
as follows: 







/ \ 






Vi 


V2 


— J 


V2 




V3 




V3 


[V4 j 







d_ 

~dt 



Following the procedure as described in the last section J can be identified as 



(27) 



J 



-1-C1W 

Here (i(t) and (2 (t) are given by 











Ut) 



1 
1 




V 



-1 + C1W oj 



(28) 



ti(t) = 2q 2 , 
C 2 {t) = -2 qi . 



(29) 



Since both qi,q 2 are determined by classical equations of motion (26), chaoticity of the 
trajectory imparts stochasticity in the dynamics of quantum fluctuations in Eq.(27). Thus, 
as elaborated in the last section, £(t) terms represent the stochastic part of the second 
derivative of the potential V(qi,q 2 ). 

If one takes into consideration of the 7-term then F(t) in Eq.(18) can be written as, 



F = F + F 1 (t) , 
14 



(30) 



where 



and 



/ 



Fi(0 



-Ci(t)rn + Ut)ri2 
The equations for the characteristic curves are 



m 

Vi 

-rji - 2<yr] 3 

-7] 2 - 27?74 






V 



(31) 



(32) 



Vl = V3 , 

m = V4 , 

m = - Ci (t)vi + (2(^2 - 27773 -vi , 

V* = C2(t)rii + &(t) 772-27774 -772 . 



(33) 



Eq.(33) describes the dynamics of quantum fluctuations in presence of dissipation. The 
mapping 77 — > rf can be found by solving the unperturbed version of Eq.(33) (i. e. the ( 
terms are omitted) for discrete small steps of r (which is consistent with the requirement 
that the correlation time is short and finite) and is given by [we refer to van Kampen^ 26 ) for 
details of treatment of multiplicative stochastic noise in Eq.(33)] 



-TVs + Vi , 



V2 T = ~ T Vi + V2 , 



27 

-r = ^ (e 2 7 r_ 1)+J74e -2 7 r 
27 



The Jacobian determinant of the transformation reads as 



(34) 
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dr) 1 
dr) 



~ e 



47T 



(35) 



where the terms of the order of r 2 are neglected. This is well within the error bound 
(a 2 r c ) as shown by van Kampen^ 26 ) and does not incorporate additional error in the analysis. 



Also note that 



dr) 



dr) T 



-47T 



(36) 



Making use of the mapping transformations 77 — > rf (Eq.34) one calculates the first and 
second derivative terms in Eq.(24) [For details we refer to Ref(26)]. The master equation 
(24) for the Henon-Heiles system can then be written down. This is 



where 



dt 



d d 

-r)3x h iVi + 2 ^3 - (C.3^1 - B 3 r) 2 + A 3 r) 2 

ar)i dr) 2 







- B' 3 rn + A 3% + B' 3 r) 2 + B 3Vl + C 3 r) 2 ) - {{( 2 (t))r} 2 - (Ci(*)>»7i)} o~ 

or) 3 

+ {r) 2 + 27774 - (((2)771 + ((1)112) + (B3V2 - C 3 ?7i - A 3 r) 2 - B' 3 r)i 

- A3771 -#s772 + ff 3 77i -C 3 77 2 )}— + 4 7 + #i 

01)4 orj^rjx 



' dr) 4 dr) 2 



+ F 1 
+ H 



d 2 



dr) 3 dr) 2 
, d 2 



+ F 2 



d 2 



dri^dr)] 



d 2 



d 2 



dr) 2 3 



dr) 3 dr), 



+ 



d 2 



dr) 4 r) 3 drjl 



(0) , 



E 1 = A A r) 2 2 + C477 2 - 77^2 (B' 4 + B 4 ) , 
E 2 = A 4 r)l + r) X r) 2 {B 4 + B' 4 ) + C 4 r)\ , 
Fi = ViV2(A 4 - C 4 ) + B 4 r) 2 2 - B' 4 r)l , 
F 2 = ViV2{A 4 - C 4 ) - B 4 r)\ + B' 4 r)\ , 

G = r) 2 2 (A 2 + B 2 ) + 77 2 (C 2 - B' 2 ) + r) x r) 2 {A 2 - C 2 - B' 2 - B 2 ) 

+ViVa{B 3 + C 3 ) + 772773(^3 - A 3 ) + 77x773(^3 - C 3 ) - r) 2 r) 4 {A 3 + B 3 ) , 

H = r) 2 (B 2 + C 2 ) + 77 2 (A 2 - B 2 ) + 771 772(^2 - C 2 + B 2 + B' 2 ) 
-rjivMs + B 3 ) + 77x773(^3 - A 3 ) - r) 2 r) 4 (B' 3 - C 3 ) 

+mvs(c 3 -B' 3 ) , 



(37) 



(38) 
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and 

A 2 = ^{{Ut)Ut-r)))e- 2 ^dT , B> 2 = J °°((Ci(t)( 2 (t - r)))e~^dr 



A, = 



fo°°((Ut)( 2 (t-T)}}e^WdT , B' s = ^{{Ciit)Ut-r)))e^rdr 



A, = 



^{{Ut)Ut-r)))rdr , B> 4 



/o° ((Ci(t)C 2 (t -r)))rdr 



> 



(39) 



B 2 = Jo 00 ((C2(i)Ci(i-r)))e- 2 ^r , C 2 



/o °((CiWCi(t-r)))e- 2 ^r 



B 3 = /o 00 ((C 2 (t)Ci(t-r)))e- 2 ^r , C 3 = / °°((Ci(OCi(* - r))) e -^rdr 



B 4 



/o°°«C2(0Ci(*-r)»rdr 



C 4 



J?°«Ci(*)Ci(«-T)»rdr 



The above equation (37) is a Fokker-Planck equation for probability distribution of quan- 
tum fluctuations for the dissipative Henon-Heiles system. It is evident that stochastic av- 
eraging over classical chaos leads to the average equation and the correlation functions 
contained in A 2 , C4 . The correlation of fluctuations of curvature of the classical poten- 
tial thus determines the drift and diffusion terms of the Fokker-Planck equation. It must 
be emphasized that this fluctuation has nothing to do with the stochasticity inherent in the 
system-heat bath model governed by Hamiltonian(2). We also point out that since the very 
notion of stochastic process in describing the curvature of the potential results in the diffu- 
sion terms, the stochastification imparts a kind of irreversibility in the evolution governed 
by the Fokker-Planck equation(37). The origin of this irreversibility is classical chaos and 
not due to any external influence. This is characteristic of the nonlinear system, itself. 



B. The solution of Fokker-Planck equation 



The appearance of the variables 771,772,773 and 774 in the diffusion terms precludes the 
possibility of an exact solution of Eq.(37). One thus takes resort to weak noise approximation 
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(this is consistent with assumption that fluctuations are not too large) scheme. The diffusion 
terms (given in the Appendix) are thus assumed to be constant. 

The resulting Fokker-Planck equation can be transformed to the following simple form 

where 

u = arii + br] 2 + ct? 3 + T] 4 , (41) 

and the constants A, A, a, b and c are given in the appendix. 

We then search for the Green's function or conditional probability solution for the system 
at u at time t given that it had the value v! at t — 0. The initial condition which is required 
to bring forth quantum-classical correspondence is represented by 

p(u,t = 0) = -e-< u - u ^ . (42) 

7T 

This means that e should be chosen in such a way that corresponds to the minimum 
uncertainty product of the initial wave packet. For notational convenience we have used 

((j>(u,t))=p(u,t) . (43) 

We now look for a solution of the equation (40) of the form 

p(u,t)\u',0) = e G(t) , (44) 

where 

G(f) = -^(u-n(t)) 2 + M*) . (45) 

We are to see that, we can, by suitable choice of Q(t), T(t) and v(t), solve Eq.(40) subject 
to the initial condition 

p(u,0)|u',0) = -e~< u - u "> 2 . (46) 

7T 
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Comparison of this with (44) with G(0) shows that 

T(0) = - , n(0) = v! , i/(0) = - . (47) 
e 7r 

If we put (44) in (40) and equate the coefficients of equal powers of u we obtain after 
some algebra the following set of equations 

r 2 dt r + r 2 ' { } 

- = -Ml , (49) 

and 

1 dv D] , . 

— = 47 - — i , 50 
udt 2T ' v ; 

where 

y = 2A , 

Di = 4A . (51) 

The relevant solution of T(t) for the present problem which satisfies the initial conditions 
above is given by 

T(t) = r(0)e- y ' + ^(l-e-^) . (52) 

7 

It is important to note that the expansion of the wave packet is determined by T(t) which 
is controlled by the two parameters, D\ and 7' which by the virtue of Eqs.(51) and (40) 
can be identified as the "renormalized" diffusion and drift coefficients, respectively. The 
origin of this "renormalization" is essentially classical chaos since these coefficients are the 
complicated functions of the correlation function of the fluctuations of the curvature of the 
potential. 
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C. Results: quantum fluctuations, expansion of phase space and entropy 



Having obtained p(r)i, r] 2 , 773, 774) we are now in a position to determine the various theoret- 
ical quantities. We calculate the quantum fluctuations of position and momentum variables. 
Since the conditional probability p is given is given by Eq.(44), this together with (48-50) 
may be employed to calculate first and second moments. Thus we express 

III I-ooP(Vi, V2, m, Va, t\v[, ^ Vs, Va, 0)mdVidV2dvsdV4 



^ III I-ooP(Vi, V2, V3, V±, tWi, rf 2 , vWa, 0)dr) 1 dr) 2 dr] 3 dr] 4 
in terms of conditional probability p. Explicit calculation yields 

Q(t) 



(53) 



Similarly we obtain 



(54) 



(55) 



The conjugate variable to 771 is r/ 3 whose average is given by 



(V3 



(56) 



Similarly 



, 1 „, N nay 
ivl) = 5? r (t) + -LL 



(57) 



Therefore the uncertainty in coordinate Ar/i and that in its conjugate momentum A 773 
are obtained as follows; 



At/ 2 = <r/ 2 > ~ (Vi) 2 = - 2 



T{t) 



At? 2 = (Vl) - (Vs? 



T(t) 



(58) 
(59) 



where the relations (54-57) have been used. The uncertainty product A771A773 at any 
time is given by 
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A Vl A m = —T(t) , (60) 

where T(t) is determined by Eq.(52) subject to initial conditions (47). This implies that we 
are to choose e = ^ to satisfy the minimum uncertainty product condition for t=0, for the 
wave packet [i.e., A771A773 = \ ] ■ 

We now discuss the following results ; 

(i) The relation (60) illustrates the evolution of quantum fluctuation as a function of 
time in terms of T(t) which by the virtue of Eq.(52) is determined by the initial condition 
T(O) [Eq.(47)] and the other two parameters D x and 7'. The early expansion of quantum 
fluctuations has been recognized as a typical signature of classical chaos on a generic quan- 
tum phenomenon^ 24 ' 25 ) . Note that Di[= AA, see Eq.(51)] is the diffusion coefficient that 
appeared in the Fokker Planck Eq.(40) [this is not to be confused with the thermal diffusion 
coefficient D in Eq.(4) which arises due to the interaction with the surroundings] and 7' 
refers to the modified dissipation rate of the system in contact with the surroundings and is 
related to A [by Eq.(51)] which is determined by Eq.(Al3). The diffusion coefficient D\ and 
modification of dissipation rate are due to the correlation of fluctuations of the curvature of 

the classical potential (±(t) and ( 2 (t) through A 2 C4 in Eqs.(39) and (Al). The origin of 

diffusion coefficient Di and the modification of 7 thus have purely deterministic origin. 

To analyze the growth of quantum fluctuations quantitatively [Eq.60] we first consider 
the dissipative classical chaotic motion governed by Eq.(26). We choose the initial condi- 
tions for energies | and |. These energy values are wellknown in the context of classical 
Henon-Heiles Hamiltonian. It is important to note that even within the restricted domain 
of weak dissipation, the dissipative Henon-Heiles system approaches a manifold reduced di- 
mensionality. The timescale over which this reduction takes place is determined essentially 
by the magnitude of the damping constant 7. This classical behaviour is illustrared in Figs. 
1 (a) and 2 (a) for 7 = 0.001 for the energies | and |, respectively. Figs. 1 (b) and 2 (6) 
depict the corresponding Poincare maps for the conservative (7 = 0.0) Henon-Heiles sys- 
tem. It is thus apparent that even weak dissipation profoundly alters the characteristics of 
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the stochastic process represented by the classical Hamiltonian chaos. The attractor clearly 
lies at the center. To calculate classical ensemble average of the quantities like (((t)) and 
{{((t)((t — t))), we carry out averaging over long time series for the given initial condition. 
The numerical procedure has been discussed earlier in Ref^ 24 ). 

Following Eq.(60) we plot the variation of uncertainty product [A^A^] (A??! and Arj 3 
are the quantum variances corresponding to position and momentum for one degree of 
freedom, respectively) as a function of time for different values of energy corresponding to 
chaotic trajectories (damping rate 7 = 0.001) in Figs 3 and 4. It has already been pointed 
out that the major input for the theoretical quantity are the chaotic diffusion coefficient D\ 

and 7' which are further related to A 2 C 4 , i.e., to classical correlation functions of the 

curvature of the potential. The theoretical curves are denoted in Figs 3 and 4 by the dotted 
lines. 

(ii) The relation (60) also shows that there exist a critical limit to the expansion of phase 
space. This limit is given by 

ATftATfel^oo = 7 r^— l ■ (61) 
2\a\cy 

The existence of this critical width is a consequence of the competetion between chaotic 
diffusion, which attempts to expand the wave packet and dissipation 7 which has the opposite 
tendency and this interplay ultimately leads to a compromised steady state. 

At this juncture it is necessary to clarify the concept of steady state of the quantum 
dissipative system as applied here. The Henon-Heiles Hamiltonian is strongly nonlinear 
and so the notion of the thermodynamic equlibrium is inappropriate. Instead we mean a 
stationary state of the quantum system in the following sense. Here we are concerned with 
an asymptotic distribution of quantum noise variables rji [zi = Zi(t) + h^rji, Zi(t) being the 
classical position or momentum variable] in terms of the probability distribution function 
(<f>(u,t)) [=p(u, t), u being the combination of quantum noise variables rji (see Eq. 41)], 
rather than a distribution of Zi(t)-s. Note that this function does not involve any classical 
contribution Zi(t) directly. The classical chaotic fluctuation in Zi(t)-s contribute to rji via 
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their classical correlation functions (in a, b, c of Eq. 41). These correlations primarily 
determine the dynamics of r\i [see (49) and (54)] through A where the role of 7 is not the 
dominant one. This implies that although the classical motion of Zi(t) settles down on an 
attractor approximately on a scale of, say, 7 _1 as noted in Figs.l (a) and 2 (a), the quantum 
noise r^-s approach the steady state very rapidly, i. e., in a few time units. Thus the 
present quantum steady state does not correspond to a settling down of classical motion on 
a attractor. Had we consider the asymptotic distribution of c-number variables Zi through 
the Wigner density function W(zi, t) , then that would have correspond to an ideal quantum 
stationary state. 

(hi) To make our analysis of irreversible evolution in presence of classical chaotic diffusion 
more quantitative, it is useful to calculate the entropy S of the Gaussian state by defining 
it (2) as 

S = -p(t)\np(t) (62) 

where p is as defined by Eq.44. In Fig 5 we show the evolution of entropy due to quantum 
noise corresponding to the classical trajectories of the dissipative Henon-Heiles system for 
the energies E — | and | and 7 = 0.001. At a very early stage the entropy change remains 
very small. It is then followed by a sharp increase and then finally tends to increase at a 
very slow rate. It is interesting to note that Zurek and Paz^ 18 ) advocated the efficacy of 
studying the evolution of entropy as a conseqence of interplay between Liouville dynamics 
and high temperature surrounding to examine the hall-mark of a nonintegrable system. 
Similar attempts had been made by us* 11 ) earlier using Husimi distribution function to 
identify the different stages of quantum evolution. 

Since the theory of stochastic fluctuations of the curvature of the potential rests on van 
Kampen's expansion in ar c as emphasized earlier, care should be taken to calculate the 
integrals (39) over the correlation functions. To implement this numerically one considers 
the first fall of the correlation functions to adjust the cut off in time for numerical evalution 
of the integrals. This is a crucial requirement for the theory which should be appropriately 
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taken care of in numerical calculation. 

We now point out a pertinent issue in the context of the nature of classical stochastic 
process considered here. It is wellknown that the Henon-Heiles model (without dissipation) 
is a typical KAM problem, i. e., it represents soft chaos and in principle never reaches the 
limit of fully developed hard chaos. This implies that, in a sense, classical noise has, in 
general, not very short correlation time. However the systematic procedure to deal with 
long correlation time when the nature of noise is rather unknown is relatively scarce. In 
principle, van Kampen's strategy as adopted here is applicable for consideration of higher 
order non-Markovian contributions. But such an extension is rather complicated both from 
analytical and numerical point of view. We therefore confine ourselves to the lowest order 
non-Markovian contribution to noise arising out of classical stochasticity and point out that 
such a description is not inappropriate in view of the short timescale over which the quantum 
fluctuations 77$ (t) determined by the correlations of the classical noise, persist. 



IV. NUMERICAL SIMULATION OF THE QUANTUM OPERATOR MASTER 

EQUATION 



For a full quantum-mechanical calculation to verify the basic theoretical propositions of 
semiclassical dynamics, we now return to Eq.(3). To solve the Eq.(3) for the Henon-Heiles 
system we choose two sets of basis vectors {|ni)} and {|n 2 )} of two different harmonic 
oscillators which satisfy (p\ 2 /2m\ + (I / 2)micu 2 qi 2 )\ni) = [(ni + l/2)huji\\ni) and (p2 2 /2m2 + 
(l/2)m 2 o;2?2 2 )|^2) = [(^2 + 1/2)^2] ^2) • The frequencies uji and uj 2 are arbitrarily adjusted 
to economize the size of the basis set. For the present purpose we choose u\ = 6.25, 
U2 = 6.20, h — 1, and 35 basis vectors. 

Quantum-classical correspondence is maintained through the construction of minimum 
uncertainty wave packets \cx quPi ) of Gaussian form in position and momentum representations 
having position (<&) and average momentum (p,) such that 
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where, 

oii = \/miVi/2[(qi) + (i / m^) (pi)] ,i = l,2 . 

The quantum evolution is followed by locating the average positions and average mo- 
menta of the initial wave packets corresponding to the initial positions and momenta of two 
classically chaotic trajectories for the two energy values | and |. Another important check 
for the numerical calculation is to keep the trace of the density matrix [Eq.(3)] equal to unity 
for the entire evolution. We have also checked that the result is robust against the variation 
of the size of the basis set. The numerical curves (solid lines) have been superimposed in 
Figs 3 and 4 for the corresponding values of energy. It may be observed that the agreement 
between the theoretical and numerical curves is quite satisfactory. This justifies the validity 
of our semiclassical approach. 



V. DISCUSSION ON THE APPROXIMATIONS, SUMMARY AND 

CONCLUSIONS 



Based on a traditional scheme of system-reservoir model we have developed a theory of 
dissipative chaotic system. We make use of appropriate ^-scaling analogous to van Kampen's 
^-expansion, of equation for Wigner quasi-probability distribution functions which takes 
into account of thermal diffusion and dissipation due to the reservoir. We have shown 
that the semiclassical approximation leads us to an equation of motion for Wigner function 
for quantum noise which is governed by the dissipation due to reservoir and the second 
derivative of the classical potential, latter being a key-point in determining the stability 
of classical motion. Since chaoticity originates from the exponential loss of correlation of 
initially nearby trajectories this derivative behaves as a stochastic (deterministic) process. 
This stochastic process is amenable to a theoretical analysis (without imposing any a priori 
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assumption about its nature) in terms of a treatment of stochastic differential equation with 
multiplicative noise. The resulting Fokker-Planck equation carries the information about 
the drift and diffusion coefficients which are expressible in terms of correlation functions of 
fluctuations of the curvature of the classical potential. As a prototypical example we have 
illustrated our analysis with the help of the Henon-Heiles Hamiltonian. 

We now make a few remarks on the approximations involved in the present treatment 
and their validity. 

(i) It must be emphasized that since the system-reservoir dissipative dynamics as gov- 
erned by the operator master equation (3) is based on Born-Markov approiximation (the 
correlation time of the reservoir must be very short (Markov) for the interaction between 
the system and the reservoir to be sufficiently small (Born/weak coupling)), the underlying 
stochastic process due to the reservoir is Markovian by construction. On the other hand 
the stochasticity due to classical chaos as inherent in the fluctuations of the curvature of 
the potential is non-Markovian since we take account of short but finite correlation time 
of this fluctuations. The construction of the associated Fokker-Planck equation is based on 
a perturbative cumulant expansion in ar c , where r c is the correlation time of fluctuation 
of the curvature of the potential. The convergence of expansion as demonstrated by van 
Kampen^ 27 ) thus allows us to retain only upto second derivative terms and as such one need 
not go to third or higher order terms to describe the dynamics. While we note that there 
exist a vast body of literature in condensed matter and chemical physics dealing with finite 
time response of the reservoir which results in frequency dependence of friction coefficient 
7, these and related aspects of dissipative dynamics are outside the scope of the present 
treatment. Our approach is similar to Graham et. al^ and Milonni et. al® in this respect. 
Thus the short time regime, we believe, is mainly controlled by the curvature of the po- 
tential and the correlation of its fluctuations or in other words the short time dynamics is 
dominated by characteristic motion of the system itself. However, in the ultimate passage 
towards equilibrium the dissipation plays a prominent decisive role. 

(ii) In the master equation(4), the classical stochasticity due to chaos and the quantum 
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noise due to incoherent processes induced by the heat bath act simultaneously and influence 
one another. We have already noted that because of ^-scaling of this equation, one arrives 
at Eq.6 in which the quantum noise term D due to surrounding does not appear in the 
lowest order. Thus it is because of strict semiclassical nature of our approach which is 
consistent with the consideration of dissipative contribution of (^f)dissipative m Eq.(4) (valid 
for kT > hu, because the propagator in the relevant integral form of the density operator 
master equation of Leggett and Caldeira [5.14 of Ref(5)] has been approximated keeping 
in view of this inequality). However a simple calculation shows that the effect of this 
incoherent contribution makes its presence felt in the next order. The quantum noise due 
to surroundings becomes appreciable only at very low temperature. 

(iii) We have already pointed out that Eq.(3) because of Born approximation is valid for 
weak damping case. We take care of this limitation by choosing small values of 7 for carrying 
out numerical simulation of quantum master equation (3) and comparing our results with 
semiclassical analysis. The latter analysis based on Eq.(4) is free from Born approximation 
and is therefore valid for both weak and strong damping limits. For a comparison over the 
entire range of dissipation one needs simulation of other kinds of master equation which are 
free from weak coupling approximation. Unfortunately most of them are not well suited for 
numerical implementation. 

We thus summarise the main conclusions of this study; 

(i) The fluctuation of the second derivative of the potential due to classical chaos is 
amenable to a stochastic description when the correlation time of fluctuations is short but 
finite. 

(ii) ^-scaling identifies an early stage of quantum evolution which is dominated by chaotic 
diffusion and dissipation but not by thermal diffusion. 

(iii) The drift and diffusion terms of the Fokker-Planck equation are intrinsic character- 
istic of dynamical properties of the system since they depend crucially on the correlation 
of the fluctuations of the curvature of the potential. The dissipation is due to the coupling 
of the system to the external reservoir which causes irreversible evolution and is truely a 
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many-body effect. On the other hand the chaotic diffusion imparts a kind of irreversibility 
in the evolution which has a strict deterministic origin and is characteristic of the nonlinear 
system itself. 

(iv) The Fokker-Planck equation is reminiscent of Kramers' equation which describes 
the Brownian motion in phase space for thermally activated processes. The Fokker-Planck 
equation also assumes a generic form for two-degree-of-freedom systems, in general. 

(v) Our results show how the initial quantum noise gets amplified by classical chaotic 
diffusion and then ultimately equilibriated with the passage of time under the influence of 
dissipation. 

(vi) We establish that there exists a critical limit to the expansion of the phase space 
which is determined by chaotic diffusion and dissipation. 

Henon-Heiles system is a classic Hamiltonian that illustrates deterministic stochasticity 
in two-degree-of-freedom systems. In view of its prototypical role played in earlier as well as 
in the present investigation, we hope that the conclusions drawn here will find quantitative 
and semiquantative applicability in other cases of dissipative two-degree-of freedom systems 
at the semiclassical level of description, in general. 
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APPENDIX A: THE TRANSFORMATION OF THE FOKKER-PLANCK 

EQUATION 



The diffusion terms corresponding to (38) under weak noise-approximation are given by 



E[ = A 4 77 2 (0) + C 4 77?(0) - 77! (0)77 2 (0) (B' 4 + B 4 ) , 
E' 2 = A^KO) + 77! (0)77 2 (0) {B A + B' A ) + CV? 2 2 (0) , 
F[ = 77! (0)772 (0)(A* - C 4 ) + S 4 77 2 2 (0) - 5^(0) , 
^2 = ^i(0)r7 2 (0)(A 4 - C 4 ) - 5 4 77 2 (0) + 5^(0) , 

G' = 77 2 2 (0)(A 2 + B 2 ) + 77 2 (0)(C 2 - 5 2 ) + 77i(0)77 2 (0)(A 2 - C 2 - B' 2 - B 2 ) 
+r h (0)rj 4 (0)(B' 3 + C 3 ) + 77, (0)773 (0)(5 3 - 40 + r/i (0)773(0) (5 3 - C 3 ) 

-^(0)774(0)^3 + ^3) , 
if' = 77 2 2 (0)( J B 2 + C 2 ) + 77 2 (0)(A 2 - B 2 ) + 77i(0)77 2 (0)(A 2 - C 2 + B 2 + B' 2 ) 
-77 1 (0)77 4 (0)(A 3 + S 3 ) + 771(0)773(0X53 - A3) - 772(0)774(0X53 - C 3 ) 
+^2(0)773(0)^3-53) • (Al) 



Zeroes in 771(0), 772(0), 773(0) and 774(0) refer to their initial values corresponding to the 
initial preparation of the coherent wave packet which is centered around the classical position 
and momentum for the chaotic trajectory. 

The Fokker-Planck equation (37) can then be written in a more compact form as follows; 




+ (771 A; + 772/ + 27773) 




+ (771 1 + 77 2 m + 27774) — (0) + 4 7 (0) + E[ 



2 d7] 4 dr] 2 




+ H\—— + — i )]<</>) 
dr) 4 dr) 3 drji 



(A2) 



where 



k = 1 - C 3 + B' 3 - B 3 - A 3 + (Ci) 
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m=l-(d) + B 3 -B' 3 -A 3 -C 3 , 
l = -B' 3 -A 3 -C 3 + B 3 -{{ 2 ) . (A3) 

We now make the following transformation; 

u — arj 1 + br] 2 + cr] 3 + r] A (A4) 

where a, b and c are constants to be determined. Using this transformation we can write 
the Eq(A2) as 

where 

A = E[ac + E' 2 b + F[bc + F' 2 a + G'(c 2 + c) + H'(c + 1) (A6) 

and 

\u = —rj 3 a - 7746 + r]ick + r] 2 c/ + 27^ 3 c + 771Z + r] 2 ^ + 27^4 . (A7) 

Making use of equation (A4) in (A7) we obtain, 

\a = ck + I , 
\b = cl + m , 
Ac = —a + 27c , 

A = 2 7 -b . (A8) 

The above relations can be used to obtain the following algebraic equation for A 

A 4 + rfiA 3 + d 2 \ 2 + d 3 \ + d = , (A9) 

where 

di = -47 , 

d 2 = 47 2 + k + m , 

d 3 = -27(fc + m) , 

d = mk-l 2 . (A10) 
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We now seek for a perturbative solution of Eq.(A9). To this end let us first note that 
in the limit 7 ^ Eq.(A9) reduces to a biqudratic form (since d\ and d 3 vanishes) whose 
solution is given by 



A = ± 



-d' 2 ± Jd' 2 - Ad 



(AH) 



where 



d' 2 = k + m . 

The lowest order perturbative solution of the the Eq. (A9) is therefore given by 

diA 2 + d 2 + 4 7 2 A 



(A12) 



A = A 



1 - 



(A13) 



4A 4 + + 2rf 2 A + d 3 _ 

For the present problem the positive real root of A is allowed which satisfies the physical 
condition (the probability distribution function must vanish at ±00). Now the values of the 
constants a, b and c which are used in Eq.(A4) can be calculated in terms of A which is given 
by Eq.(A13) as follows 

b = 2 7 -A , 
I 



A(2 7 -A)-A; ' 
a = (27 - A)c . 



(A14) 
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Figure Captions 



1. (a) Plot of qi vs pi on the Poincare suface of section for the Henon-Heiles sytem with 
damping constant 7 = 0.001 and initial energy E — |. (b) Same as in (a) but for 
7 = 0.0. 

2. (a) Same as in Fig. 1 (a) but for E — |. (6) Same as in (a) but for 7 = 0.0. 

3. Plot of uncertainty product (A771A773) with time for the system as in Fig.(l). The 
continuous line represents the numerical simulation of the master equation (fully quan- 
tum). The dotted line refers to semiclassical calculation ( Eq.(60)). ( Both units are 
arbitrary) . 

4. Same as in Fig. (3) but for E — |. 

5. Plot of evolution of entropy with time for damping constant 7 = 0.001. (a) E — | and 
(b) E=\. 
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